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Abstract - It is shown using dimensional analysis that the maximum current density Jqcl 
transported on application of a voltage V g across a gap of size D follows the relation Jqcl ~ 
h 3 ~ 2a V£" J D 5 ~ 2a . The classical Child-Langmuir result is recovered at a — 3/2 on demanding that 
the scaling law be independent of h. For a nanogap in the deep quantum regime, additional inputs 
in the form of appropriate boundary conditions and the behaviour of the exchange-correlation 
potential show that a — 5/14. This is verified numerically for several nanogaps. It is also argued 
that in this regime, the limiting mechanism is quantum reflection from a downhill potential due 
to a sharp change in slope seen by the electron on emerging through the barrier. 
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Introduction. — The Child-Langmuir (CL) law [T] 
is the cornerstone of classical space charge limited flows. 
It determines the maximum current (Jcl) that can flow 
across a gap of size D on application of an electric held 
'—Vg/D, given the existence of mono-energetic electrons 
at the cathode [HGIj. The limiting mechanism is fairly 
straightforward: the mutual repulsion between the elec- 
trons gives rise to a barrier in the net electrostatic poten- 
tial and an energy spread in the initially mono-energetic 
beam of electrons. Beyond the limiting current, the ex- 
cess electrons have energy slightly lower than the barrier 
height and arc thus reflected back to the cathode. This 
picture holds for electrons with zero energy as well (rela- 
tive to the cathode potential) where a more common, al- 
beit equivalent, criterion of the limiting mechanism is the 
total cancellation of the applied held by the space charge 
held at the cathode and can thus be arrived at from the 
Poisson equation of electrostatics. 

The quantum treatment for studying steady state 
charge transport across a gap is based on a self-consistent 
solution of the Schrodingcr and Poisson equations [H[S]. 
As in the classical case, there exists a maximum current 
in this formalism beyond which a solution to the coupled 
set of equations ceases to exist. The exact limiting mech- 
anism is however largely uninvestigated with the studies 
largely centred around scaling laws. 

The purpose of this letter is to (i) provide a unihed scal- 
ing law across all regimes whenever a power law behaviour 
holds, (ii) probe the deep quantum regime using additional 



inputs such as the boundary conditions and the behaviour 
of the exchange-correlation potential and (iii) investigate 
the limiting mechanism in the quantum domain. 

The model that we study pre-supposes the existence 
of free mono-energetic electrons at the cathode, whether 
due to thermionic, held or photo emission. As in case 
of the classical Child-Langmuir law, the quantum frame- 
work merely seeks to predict the maximum steady current 
that can be transported on application of a field. Impor- 
tantly, the classical or quantum Child-Langmuir law does 
not in itself address the question of emission. Together 
with emission laws however, the limiting Child-Langmuir 
current plays a role in understanding current- voltage char- 
acteristics. 

Quantum effects become important at low applied volt- 
ages or small gap size and give rise to voltage scaling other 
than the classical three- halves law (V g ), a subject of 
considerable interest due to the current focus on nanos- 
tructured materials. In the deep quantum limit, the small 
phase volume permits few electrons so that the Hartrce 
potential is inconsequential and the exchange correlation 
potential dominates. We shall focus on this regime in this 
communication and try to understand the limiting mech- 
anism and explore the scaling behaviour with respect to 
the applied potential. 

The quantum phenomenon of reflection, tunneling and 
the fermionic nature of electrons are best incorporated 
within the Kohn-Sham density functional theory [B] with 
an effective potential, V e ff, that includes a parametrized 
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form of the exchange-correlation potential (such as due to 
Pcrdcw and Zunger [7J) obtained within the local density 
approximation (LDA). Since, the emission mechanism is 
not considered here, limitations of LDA, such as the fact 
that it does not correctly reproduce the image potential, 
is not of much consequence. 

The basic inputs for implementing Kohn-Sham theory 
to study transport across nanogaps are twofold. The first 
concerns the boundary conditions for the Hartree poten- 
tial at the gap boundaries. Since the applied voltage differ- 
ence is proportional to the difference in chemical potential 
between the injection and collector planes, the exchange- 
correlation potential at the two ends play a significant role 
in determining the boundary conditions in the deep quan- 
tum regime and its neglect can lead to erroneous results 
[5]. The second input lies in the determination of the 
boundary conditions for the wavefunction and is related 
to assumptions for the potential beyond the gap for pur- 
poses of matching the wavefunction. The simplest of these 
is to assume that the effective potential assumes a constant 
value (abruptly) on either side of the gap. This approach 
is however known to give misleading results when the 
exchange-correlation potential is neglected altogether [5] . 
First, an abrupt change in potential for the Schrodinger 
equation can result in considerable quantum reflection at 
low applied voltages or small gap size. Further, the cou- 
pling of the Schrodinger and Poisson equations makes the 
boundary conditions for the wavefunction and Hartree po- 
tential inconsistent [9] . The problem can be addressed by 
simply assuming that the transported charged particles 
are not reflected back into the gap and that the effective 
potential varies smoothly across the interface. An approx- 
imate transmitted wavefunction is then readily available 
within the semiclassical approximation, which can then 
be matched with the wavefunction in the nanogap at the 
collection plane. 

Formalism. — We shall first familiarize ourselves with 
the equations to be solved self-consistently and subse- 
quently arrive at the boundary conditions that they need 
to satisfy. The effective potential energy, V e /f, in the 
Schrodinger equation 



d 2 ifj/dx 2 + V e ffip = EiP 



(1) 



takes the form V e fj = —eV + V xc x Eh, where Eh = 
e 2 /(47reofflo) is the Hartree energy and ao the Bohr radius. 
The exchange correlation potential within LDA takes the 
form V xc = e xc - (r s /3)de xc /dr s , where r s = [3/(47m)] 1/3 
is the Wigner-Seitz radius, n the electron number density 
and e the magnitude of the electronic charge. In the above, 
V is the Hartree potential satisfying the Poisson equation 



d 2 V/dx 2 = en(x)/e = e\tp(x)\ 2 /e 



(2) 



as e xc = e x + e c where e x = — (3/27r) 2 / 3 (3/4r s ) is the ex- 
change contribution for a uniform electron gas while e c , the 
correlation contribution, is represented by a parametrized 
form of the random phase approximation result. Un- 
less otherwise specified, we shall restrict ourselves to the 
Pcrdcw-Zungcr [7J paramctrization. 

We first note that the applied voltage difference, V g = 
— {HC — ^e)/b where fie an d (J-e refer respectively to the 
chemical potential at the collector and injection planes. 
For convenience, we consider the reference as [Ie = 
-eV(0) + V xc {0) x E H = so that E = refers to in- 
jection from the Fermi level. Thus V(0) = V xc (0) x Eh fe- 
lt follows that the chemical potential at the collector 
is -eV(D) + V XC {D) x E H = -eV g . Thus V(D) = 
V g + V XC (D) x Eff/e. In writing the above, we have im- 
plicitly assumed continuity of the chemical potential at 
the interfaces under steady-state conditions. 

The coupled Schrodingcr-Poisson system with the 
above boundary conditions for V are to be solved self- 
consistently under the assumption that a current density 
J is transported across the collection plane. This accounts 
for one of the two boundary conditions required to solve 
the Schrodinger equation. The other is derived under an 
approximation. 

The simplest of these is the assumption that the effec- 
tive potential is constant beyond the collection plane. As 
discussed in the introduction however, this can lead to er- 
roneous results, mainly because a self-consistent solution 
requires Vf*JD) = 0. If the exchange-correlation contri- 
bution is small, this implies V'(D) = 0, a condition that 
is found to be violated except in regimes where a classical 
description provides an adequate solution. 

At the next level, assuming that there is no reflection 
back into the nanogap and that the effective potential 
varies sufficiently smoothly at the interface, the trans- 
mitted solution can be expressed as a first order WKB 
wavefunction 
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>£ r p(x')dx' 



VP(x) 



(3) 



where p(x) = ^2m(E — V e ff(x)) is the classical momen- 
tum. The constant C can be fixed in terms of the known 
quantities J, E and V e ff(D). The approximate transmit- 
ted wavefunction, iptrans can then be used to match the 
gap wavefunction ?p and its derivative at the boundary. 

The coupled system of equations can be solved by first 
expressing the complex wavefunction 



1& = ^T q(x)e W ^ 



(4) 



with boundary conditions that we shall shortly specify. 
The exchange correlation energy density, e xc is expressed 



in terms of a real amplitude q(x), phase 6{x) and the 
characteristic density no = 2eoV g /3eD 2 . Using dimen- 
sionless normalized variables x = x/D, (f> = V/V g , A = 
D/Xo where the electron de Broglie wavelength Ao = 
y / h?/2meV g , e = E/eV g , J = J/Jcl the Schrodinger and 
Poisson equations can be expressed respectively as [S] 
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where <j) xc = 
(r s /3)de xc /dr s 



A J 2 
h-c - t— t]q = and (5) 



9q 4 



dx 2 3 y 

V XC /(j)g, (f>g = eVg/E H , V x , 



(6) 



and 



dV xc 
dx 



it follows that 



dV XC l 
dr. 



dr s 
|..(i) x — |x=l 
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(7) 
(8) 



with r s = [3A/20 9 g(x)] 2 / 3 , 7 = -0.1423, ft = 1.0529 and 
ft = 0.3334. The parametrized form for the correlation 
energy density in Eq. [S] is due to Perdew and Zunger [7] 
with the parameters obtained by fitting to the random 
phase approximation results. 

With the above scaling, the boundary conditions for the 
scaled Hartree potential (Eq. [5]) are 



0(0) 

0(1) 



0zc(O) 

1 + 0.T C (1). 



(9) 
(10) 



The boundary conditions for q(x) can be derived under 
the assumption that the transmitted wavefunction is given 
by Eq. [3] and carries a current density J. Using Eq. [3] and 
J = eih(ip*ip' - V>"0*')/ 2 ™: ii; follows that \C\ 2 = Jm/e. 
Without any loss of generality, we shall assume that C 
is real. Thus C = y/Jm/e. On equating ifj (Eq. E]) and 
iptrans (Eq. [3J) at x — D (x — 1) , we have 



<?(1) = 



(l) + e-<Ml)] 1/4 



11/4 



(11) 

(12) 



<?'(!) = — 



where 



C = i + 



(1)1 

[1 + e] 5 / 4 c 



-(3X/2) 2 / 3 (dV xc /dr s )\ rs{1) 



(16) 



(17) 



[1 + e] 5 / 4 

As V xc — > 0, C — > 1 so that q'(l) takes the value predicted 
on ignoring exchange-correlation effects altogether [5] . 

Voltage Scaling: Dimensional Analysis. As the 

Schrddingcr equation dictates the deep quantum regime, 
the scaling with applied voltage should differ from the 
Poisson equation induced V g scaling. The existence of a 
power law behaviour in the quantum regime cannot be as- 
sumed outright. However, numerical results clearly show 
a power law behaviour in the deep quantum regime (see 
Fig. []}. 




where Eq. [T^] follows on using the boundary condition 
0(1) = 1 + <f> xc (l). In arriving at Eq. [T2J we have equated 
9(1) = (1/fr) J p(x')dx', in effect fixing the value of the 
reference point x r . 

The boundary condition for dq/dx can be similarly de- 
rived by equating the derivatives of ip and iptrans- How- 
ever, since <j) xc depends on q(x), the value of q'(l) must be 
computed from the equation 



</(!) = -7 



1 /2J [0'(l)-V x ' c (l)/0 g ] 



1 



<?Mi)] 5 / 4 



(13) 



where the primes indicate differentiation with respect to 
x. Since 



Fig. 1: The maximum scaled current density J vs V g for three 
different nanogaps (i) D = lOnm - top set of points marked 
as diamonds (ii) D = 30nm - solid squares in the middle (iii) 
D = 50nm - the bottom set marked by solid circles. The best 
fitting straight line in the deep quantum regime is also shown. 
In all three cases, the slope is in the range -1.16 to -1.14. Note 
that the voltage at which the deep quantum regime sets in is 
smaller for larger gaps. 

Expressing the scaling in terms of dimensionless quan- 
tities as J ~ (V g /V s )P where J = J/Jcl and the voltage 
scale V s = h 2 /(2meD 2 ), it follows that 



epeP+^h- 2 ? F/ +3/2 

m l/2-/3 £)2-2/3 ■ 



(18) 
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Writing j3 + 3/2 = a, it follows that 



J~ e - 



h 3 ~ 



VI 



D 5 - 



(19) 



giving the scaling J ~ V" / D 5 ~ 2a as shown in [9]. 

The scaling law above is based on numerical observa- 
tions. It can however be derived from purely dimensional 
analysis. The relevant equations are the Schrodinger 
(Eq|T]) and Poisson equations (Eq. [2]) together with the 
current equation 



2im 



(20) 



Redefining ijj = ^>/v/eo, it follows that the Schrodinger and 
Poisson equations are free of eo while the current density 
takes the form 



J 



Aim 



i>*^i>]. 



(21) 



Thus, i\) no longer depends on eo so that J ^ eo. However 
-0 continues to depend on e, h and m and hence their 
dependence on J need not be linear. We are now in a 
position to implement standard dimensional analysis to 
ascertain the scaling law. 

Since, the critical current can in general depend on the 
quantities V g , D, eo, fi, m and e, we demand 



j-^V° i D 6l H 62 e m 53 e 6i ~ 



= C 



(22) 



where Co is a dimensionless constant. On expressing all 
the quantities in terms of mass, length, time and current 
(M, L, T, A), and equating powers of these to zero, we ob- 
tain four equations 



5± — a 


= -1 


(23) 


2S 2 + 2a + S 1 = 


= 1 


(24) 


64 — 62 — 3a = 


= -4 


(25) 


$2 + £3 + a - 


= 1 


(26) 



in terms of five unknown quantities. Expressed in terms of 
a, the exponent of V g , we obtain 61 = 2a - 5, 82 = 3 — 2a, 
Ss = a — 2 and 64 = a — 1. Thus, standard dimensional 
analysis gives us Eq. [TO] confirming that any power law 
behaviour must follow this relation. 

Note that the Child-Langmuir law follows automati- 
cally on demanding that the exponent of H be zero. Thus 
a = 3/2 so that J - V g ,2 /D 2 [Mil]. 



Scaling in the deep quantum regime. While 

dimensional analysis leads us to the general scaling law 
j _ fp-2 a y<x J D $-2 a ^ the value of a b se rved in Fig. [T] 

can be understood by the following analysis that involves 
additional inputs. We first note from Eq. (fi"2"j) that 



J 



? 2 (1)- 



To determine the dependence of g 2 (l) on V g , we shall in- 
vestigate the quantity 



dV, 



<z 2 (i) 



dV 



q 2 {x)\x= 



(28) 



where the equality follows on using the alternate but 
equivalent boundary conditions for the electrostatic po- 
tential: 



V(x = 1) 






(29) 



V(x = 0) = -^(T4c(0)-F xc (l)) (30) 

with the injection energy redefined equivalently as E = 
V xc (1)Eh- Note that this does not result in any change 
in the expression for q(l) or q'(l) since </>(l) + e — <j) xc (l) 
remains invariant (equals unity) for both sets of boundary 
condition and injection energy. 
Further, since 



&<*-»-&/±n~ 



we have 



dV r 



■q 2 (l) = 2q{l) q \l)/V\l) 
J 



2g(lV(l) 

v g <t>'{i) 



(31) 



(32) 



(33) 



We now note that in the quantum regime £ >> 1 so that 
C ~ \ 2/3 y/J(dV xc /dr s )\ ra{1) . As A - V g 1/2 , it follows that 



dVg 



v7 



V 9 \ dr s )\ r3 (l) 



(34) 



Note that the electron density is small at x = 1 compared 
to x = and is negligible at low applied voltages. Thus 
r s (l) is large in quantum regime. It can be verified that 
for the Perdew-Zunger exchange-correlation potential, 



dV^ 
dr„ 



k(i) 



,2/3 



(JJcl) 2/3 ~ J 2/3 V & 



(1) ~ J 2/3 (35) 



Thus, 



dV, 



<? 2 (1) 



Va /3 JV 6 



(36) 



(37) 



Assuming a power law behaviour J ~ V" , the above equa- 
tion translates as 



ya'-l+a'/6 ^ y-7/3 



(38) 



(27) On equating the powers, we have 7a'/6 = —4/3, or 
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-8/7. 

5/14 



(39) 

Thus, J ~ Vg*'' or J ~ Vg 0/1 * so that a = 5/14. The 
value of a' = —8/7 ~ —1.143 agrees very well with the 
numerical results presented in Fig. Q] where the exponent 
a' is found to lie between -1.16 and -1.14 in all three cases. 




Fig. 2: Typical effective scaled potential curves for D = 30nm 
and V g — 2 X 1CP 7 V in order of increasing current density from 
bottom to top. The horizontal dashed line is cpeff = 0. The 
values of J are 20,30,35,37 and 39.85 (topmost). The critical 
scaled current density is around 39.9. 

The Limiting Mechanism. In the classical 

regime, the Schrddingcr-Poisson system decouples and the 
Poisson equation alone dictates the limiting behaviour. 
In the deep quantum regime of low applied voltages in 
nanogaps, numerical investigations show that <j>(x) is lin- 
ear and follows the applied voltage curve. Thus, the 
Schrodingcr equation alone decides the existence of a phys- 
ically acceptable solution. We shall briefly consider the 
limiting mechanism for the model considered in this com- 
munication. 

In the steady state, an acceptable solution is a smooth, 
spatially non-oscillatory behaviour of the effective poten- 
tial. At very low applied voltages in a nanogap, there 
are too few electrons to alter the electrostatic potential 
curve from the linear applied potential significantly. Su- 
perimposed on this however is the repulsive exchange- 
correlation potential which decides the effective potential 
in the Schrodinger equation. Typically, <f> xc is small at 
large r s and falls below the electrostatic potential close to 
the right boundary. Thus the attractive potential dom- 
inates near the anode. Near the cathode however, r s is 
small and <j) xc dominates over <fi. Typically, an acceptable 
limiting effective potential in the deep quantum regime has 
a barrier towards the cathode, passes through zero and be- 
comes attractive closer to the anode. For D = 30nm and 
V g = 2 X 10~ 7 V, the change in scaled effective potential 
with current density is shown in Fig. [5] The topmost po- 
tential curve is close to the critical current density. Note 
that such a solution ceases to exist beyond criticality. 



Two observations can be made immediately from the 
change in potential curves as criticality is approached. 
First, the barrier height and width increase as the cur- 
rent density is increased. Second, in the region where the 
electron is free, from the point where it tunnels through 
the barrier to the collector plate at x = 1, the change in 
slope increases as criticality is approached (see Fig. [2]). 

The increase in barrier height and width however can- 
not be the limiting mechanism for the model considered 
since even a tiny transmission coefficient cannot violate 
any of the boundary conditions, in particular that of an 
outgoing transmitted wave carrying a current density J 
to the collector. It of course requires a large number of 
electrons to be present at x = 0, due, for example, to 
photo-emission. 
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Fig. 3: The maximum slope of the effective potential in the free 
region vs V g for three different nanogaps (i) D — lOnm - top set 
of points marked as circles (ii) D = 30nm - solid squares in the 
middle (iii) D = 50nm - the bottom set marked by diamonds. 
The best fitting straight line in the deep quantum regime is 
also shown. In all three cases, the slope is in the range 0.33- 
0.34. The maximum slope is measured in the region where the 
electron is free. In the deep quantum regime, the maximum 
slope occurs at the point where (f> e ff becomes zero while in the 
classical regime, the maximum slope shifts to x = 1. 

A sharp change in slope of the effective potential in the 
downhill region is however known to cause quantum re- 
flection [T3"irB] just as a sudden change in refractive index 
leads to optical reflection. Such a phenomenon does occur 
in the approach to criticality as seen in Fig. [5] In order 
to test the hypothesis, we study the dependence of maxi- 
mum slope on the applied voltage for different nanogaps in 
the free region. Fig. [3] shows a universal behaviour in the 
quantum regime of low applied voltages. The dependence 
of maximum slope on the applied voltage is clearly identi- 
cal in all three cases with \dcj) e f f / dx\ ~ V„ s with S ~ —0.33 
[15] . Since the slope d(f> e ff/dx — > as x — > 1, the onset of 
criticality appears to be linked to the change in slope seen 
by the electron as the electron emerges from the barrier. 
Note that a crucial input in deriving the boundary con- 
ditions is the assumption that an outgoing wave carries a 
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current density J . This breaks down in the presence of 
quantum reflection from a downhill potential and is pos- 
sibly the limiting mechanism |16j . 

Summary and Conclusions. — We have considered 
the problem of maximum transmitted current on applica- 
tion of a bias, assuming that a sea of electrons exist at 
the end having lower potential, for example due to photo- 
emission. The quantum model considered is 1-dimcnsional 
and includes exchange-correlation effects. The applied 
bias is carefully fixed in terms of the difference in chemical 
potential at the two ends and the boundary conditions for 
the wavefunction amplitude are derived using a WKB ex- 
pansion at the collector plane to account for a consistent 
solution of the Poisson-Schrodinger system. 

It was shown using standard dimensional analysis that 
the scaling law applicable across all regimes is Jqcl ~ 
h 3 - 2a V ( f J D 5 - 2a . The classical Child-Langmuir law fol- 
lows on demanding that the exponent of h be zero. Thus 
a = 3/2 and J QCL ~ V g 3/2 /D 2 . 

It was found numerically that for the Perdew-Zungcr 
exchange-correlation potential: (a) the maximum trans- 
mitted current falls below the classical prediction till the 
quantum regime is reached (b) the applied voltage sig- 
nalling the onset of the quantum regime decreases as the 
gap size is increased (c) in the deep quantum regime, the 
maximum transmitted current density is higher than the 
classical prediction and the region is also marked by a uni- 
versal behaviour in the scaling of the current density with 
applied voltage. The voltage scaling was shown analyti- 
cally to follow J ~ V" with exponent a = 5/14. The ex- 
ponent was found to be close to the numerically observed 
values. 

We conclude this communication with a few remarks: 
(i) the results presented here are directly applicable ex- 
perimentally only when electrons are made available at 
the emitter end by an alternate mechanism such as photo- 
emission. In case of field emission, these results can only 
provide an upper bound to the transmitted current density 
(ii) the deep quantum regime depends sensitively on the 
exchange-correlation potential and predictions of higher 
transmitted current and universal behaviour may vary 
with a different parametrizations of V xc . It may be hoped 
that careful experimental observations will help decide the 
appropriate form of of correlation energy applicable for 
electron transport in nanogaps. 
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